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Abstract 

~ T J^ "Lorentz Invariance Violation (LIV) may be a good observational window on Quantum Gravity physics. Within last few years, 

^—1 all major gamma-ray experiments have published results from the search for LIV with variable astrophysical sources: gamma-ray 

<^^ bursts with detectors on-board satellites and Active Galactic Nuclei with ground-based experiments. In this paper, the recent time- 

£n of-flight studies with unpolarized photons published from the space and ground based observations are reviewed. Various methods 

C\ used in the time delay searches are described, and their performance discussed. Since no significant time-lag value was found 

D within experimental precision of the measurements, the present results consist of 95% confidence level limits on the Quantum 

l"H Gravity scale on the linear and quadratic terms in the standard photon dispersion relations. 

^H ■ 

fvi Keywords: Lorentz Invariance, Quantum Gravity, Photon propagation, Gamma-ray bursts, Active galaxies 
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1. Introduction 

Postulating various symmetries in Physics is largely justi- 
fied by the quest for unified theories which rely on mathemat- 
ical symmetries and on invariance of the fundamental laws. In 
the far past, the early Universe was dominated by a dynamical 
symmetry evoluting to the complex diversity of broken symme- 
tries. In the history of the Universe the unique theory of Space 
and Time was characterized by Plank scales: /p = (G/z/c 3 ) 1 '' 2 » 
1(T 33 cm, f P * 1(T 43 s and E v « 1.2 X 10 19 GeV. Later, the 
emergence of new phenomena on hierarchical levels lead to our 
present world composed on one side of galaxies of stars and 
on small scales of atoms, molecules and life. At present there 
is no consensus on the correct approach to the unification of 



small and large-scale physics: the nature of the small-scale of 
the space-time could be derived from the idea of "fuzziness" 
(following to the Hei senberg uncerta i nty principle) or quan 



turn foa m discussed by Wheeler! ( 1982r)lAmelino-Camelia et al 
dl998al) or lAm elino-Camelia (2002), or sharply defined set of 



discrete points (JGambini and Pullini 119991) . This "lumpiness" 
or "discreteness" may affect photon propagation in the vac- 
uum resulting in propagation anomalies. In consequence, the 
speed of light would differ with energy when photons travel 
through large distances. The searches for a physical theory uni- 
formly valid in large (strong field gravity, High Energies, cos- 
mology) and small scales (atoms, nuclei, elementary particles) 
need empirical tests across these scales. The quickly develop- 
ing domain of studies with astronomical observations of distant 
sources may shed new light on the subject. 

The purpose of this paper is purely experimental and 
no critical discussion about various theoretical approaches 
will be presente d here, a s it m ay be found in e.g. 
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Amelino- Camelia and Smolinl (120091 and references therein). 
The formalism in use will be limited to the one adopted in the 
time-of-flight studies with unpolarized photons from distant as- 
trophysical sources. These constraints are usually translated 
into Quantum Gravity Scale limits in the frame of considered 
models. In consequence, the naming of Quantum Gravity will 
be used as a generic appellation of the quantum effects in the 
space-time structure. 

The constraints on Lorentz Invariance Violation (LIV) 
derived here are much less competitive than those ob- 
tained in birefringence studi es in the frame of the Standar d 
Mod el Extensions scheme of Colladav and Kosteleckvl (11998); 
Kosteleckv and Russelll (120081) . They should be considered as 
a complementary way to study the LIV effects and their conse- 
quences in the most model-independent way. 

The Lorentz Symmetry has been successfully tested with 
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experiments dealing with particle physics domain. However, 
there was no confirmation yet of its validity at large scales. 
Large theoretical interest in possible high energy violation of 
local Lorentz Invariance (LI) in the past decades was driven 
by possible hints for the Quantum Gravity (QG) as a lo- 
cal LI may not be an exact symmetry of the vacuum. As 
stated above, there are several reviews available on the subject 
dAmelino-Camelia and Smolinll2009tlEllis et alll2009h . 

The present review paper is organized as follows: the next 
section introduces the search for LIV and QG scale limits 
with astrophysical sources and briefly describes the formalism 
in use. Section [3] presents the results obtained with Pulsars, 
Gamma-ray Bursts and Active Galactic Nuclei. Finally, results 
from latest analyses are summarized and discussed in section|4] 

2. Tests of Lorentz Invariance and search for Quantum 
Gravity with Astrophysics 

2.1. Testing Lorentz invariance with astroparticles 

Several astrophysical messengers traveling from distant 
sources can be considered when testing LIV: photons, elec- 
trons, Ultra High Energy Cosmic Rays (UHECRs), neutrinos 
and gravitational waves. Each type of messenger presents its 
proper advantages and drawbacks. The propagation of pho- 
tons can be simply described with dispersion relation formu- 
las discussed below. On the other hand the detection of the 
high energy photons suffers from limitations in the energy lever 
arm due to the so called "gamma horizo n" effect related to the 
intera ctions with infra-red background (B lanch and Martinez , 
2005). The electrons emitted by pulsars synchrotron radi- 
ation allow very precise measurements of t he time delays 
but their signals come from close-by sources (Jacobsonet al., 
20061: iLiberati and Maccionel 2009; St ecker and Scully, 2009). 



The UHECRs are currently used to investigate the possible 
modifications of the Greisen-Zatsepin-Kuzmin (GZK) thresh- 
old responsible fo r a cut-off in the UHECR spectrum at ~ 
10 18 ^ eV (see e.g. iGalaverni and SiglL 120081: iMaccione et al. . 
20091) . These studies provide already very stringent limits on 
LIV but are also very dependent on the particle physics model- 
ing and on the anisotropy of their emission. The neutrinos from 
the astrophysical sources as well as the gravitational waves have 
not been yet observed. However, in case of a co incident detec- 
tion, the results on LIV w ould be very stringent (Pradier, 2009; 
van Elewvck et all 200 9J). As recently remarked by several au- 
thors (see e.g. lAmelino-Camelia and SmolinL |2009), the LIV 
parameters may vary for different messengers as protons, elec- 
trons, gamma, neutrinos and gravitons, depending on a consid- 
ered theoretical frame. 

LIV studies wi th photons from cosmol o gical so urces were 
first proposed by lAmelino-Camelia et al.l dl998al) . As sug- 
gested in this article, the tiny effects due to quantized space- 
time can add up to measurable time delays for photons from 
distant sources. The hypothesis of energy dispersion could be 
best verified in sources that show fast flux variability, are at cos- 
mological distances and are observed over a wide energy range 
between keV and TeV This is the case of Gamma Ray Bursts 



(GRBs) and Very High Energy (VHE) flares of Active Galactic 
Nuclei (AGNs). Both types of sources are the preferred targets 
of these "time-of-flight" studies which provide the least model 
dependent tests of the Lorentz Symmetry breaking. The case 
of pulsed emission by Pulsars has also been considered, and 
provide valuable results as discussed below. 

The measurement of the arrival time of photons in the de- 
tectors may only be possible if each trigger acquires a precise 
time-stamp provided by a spatial global positioning receiver. 
The precision in time below a micro-second, which is a key 
requirement of the described studies, allows a production of 
precise light curves for a given source on a few second time- 
scales. In particular, the Galileo network[j which is a global 
navigation satellite system (GNSS) currently being built by the 
European Union (EU) and the European Space Agency (ESA), 
will provide in future the needed accuracy of the detection time 
of individual photons. In addition, the UTC time provided by 
the GNSS makes possible the multi-wavelength analyses which 
combines results from various space experiments operating at 
lower energies such as SWIFT, INTEGRAL or the future satel- 
lite SVOM and higher energy satellites such as Fermi and AG- 
ILE, all integrated in the worldwide alert system (GCN). 

The Figure of Merit of a given source for the time-lag 



measurement may be formulated as dAmelino-Camelia et al 
Il998bh : 



L E 1 



1 = ~ 



c Eq G At' 



(1) 



where L is the distance of the source, At is the measured time- 
lag and Eqg is the Quantum Gravity scale. 

2.2. Measurement of the propagation effects with photons 
Following the formalism employed by Ijacob and Piran 



(2008), a general model in which the Lorentz Symmetry is bro- 
ken at a very high energy denoted in the following by Eqg, is 
assumed. A natural value for the QG energy scale would be 
the Planck energy, however a large range of values should be 
considered. 

For photons with much lower energies than Eqg, the first 
and second order corrections to the energy-momentum relation 
are generally taken into account, as shown with n = {1,2} in the 
following equation: 



E*-p 2 c^±pV\f- 



(2) 



Considering sources at the cosmological distances, the anal- 
yses of time lags as a function of redshift requires a cor- 
rection due to the expansion of the universe, which depends 
on the cosmological model. Following the analysis of the 
BATSE da ta and more recently of the HETE-2 and SWIFT 
GRB data dBolmont et allbOOalEms et alll2003ll2006ll2008l) 
and within a framework of the Standard Cosmological Model 
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Figure 1: Light curves of the 15 GRBs detected by HETE-2 in the energy range 6^400 keV used in Bolmont et al. 12008). The bursts are sorted by increasing z from 
top left to bottom right. While z increases, the signal-to-noise ratio decreases. X-axes are graduated in seconds. 



(Bah call et allll999l) with flat expanding universe and a cosmo- 
logical constant, the difference in the arrival time of two pho- 
tons with energy difference AE is: 



, AE 
for a linear term and 



z (1 + z)dz 
Hz) 



Af = 2 H o' 



AE 1 



I 



(1 + zfdz 
Hz) ' 



(3) 



(4) 



for a quadratic term, where AE 2 is the quadratic energy dif- 
ference. h(z) is given by h(z) - \/Qa+Q^(1+zP. Q m , Qa 
and Ho are parameters of the Standard Cosmological Model 
(Q,„ = 0.3, Q A = 0.7 and H = 70kms"' Mpc" 1 ). The time-lag 
may decrease or increase with AE depending if the model lead 
to sub-luminal or super-luminal photons. 

For completeness, it is necessary to mention that the for- 
malism of Equations [3] and |4] even if it has been commonly 
used by gamma-ray astrophysicists for more than a decade, 
is not the only one u sed in the community studying LIV. 
iKosteleckv and Mewesl (2008) for example write these expres- 
sions as|3 



A^ = y r z d + z)dz 

AE 'J,, U h(z) 



for a linear term and 



At 

AE 1 



Jo Ho 



zfdz 



Hz) 



(5) 



(6) 



for a quadratic term and set the limits on parameters Y\ and Y q . 
These parameters can easily be deduced from E GG and EJL, by 



Yi = 



E l 



and 



Y a = - 



2 (Eqc) 2 ' 



(7) 



(8) 



In the following, and as it is straightforward to switch from one 
formalism to the other, the one of Eqs.|3]and|4]will be used and 
the limits will be given on Eq G and El G . 

In order to probe the energy dependence of the velocity of 
light induced by LIV, analyses of time lags studied as a func- 
tion of redshift with several sources are the best way to eval- 
uate effects due to the photon propagation. In some cases, the 
parameter Ki is defined as follows: 



K,(z) 



£ 



(l+z)dz 
Hz) ' 



(9) 



to take cosmological effects into account, except for nearby 
sources for which the simple figure of merit (Eq.Q]) is a good ap- 
proximation. Till 2007 and the publications of Jacob and Piranl 



2 Refer to Eq. 145 of Kostelecky and Mewes 12008) for the exact expression 
of Y, and Y q . 



(2007, 20081) and lBolmont et al. (2008), the formula in use in all 
presented results had a term (1 + z) missing under the integral, 
thus yielding underestimated values of limits on the Quantum 
Gravity scale. 

The energy-dependent time-lags may arise from propaga- 
tion effects in the LIV scheme discussed in this paper. On the 
other hand, possible effects intrinsic to the astrophysical source 
emission could also produce energy-dependent time-lags. The 
analysis as a function of the redshift ensures in principle that the 
results are independent of source-induced time-lags. In the first 
order of the dispersion relations, the evolution of the time-lags 
as a function of z can be written as: 



(At)=aKi(z) + b(l+z), 



(10) 



where a and b parameters stand for extrinsic (Quantum Gravity) 
and intrinsic effects, respectively. When the redshift study is not 
possible, the intrinsic (or source) effects are treated in average 
or assumed to be negligible. This is the case for most results 
presented in Section[3] 

The most expanding domain of LIV tests at the moment 
is related to the astrophysical experiments with photons. The 
space missions are equipped with excellent detectors for vio- 
lent event detection such as pulsar and GRB emission. They 
included BATSE, HETE-2 and INTEGRAL in the past, SWIFT 
and Fermi at present. The ground based telescopes H.E.S.S., 
MAGIC, VERITAS and CANGAROO provide the highest AE 
with AGN flares but are limited in time variability. 

2.3. Methods 

There are various methods for the time-lag determination 
and subsequent derivation of the Quantum Gravity scale. In or- 
der to measure a tiny deviation from its standard value of the 
light velocity and reach 10 19 GeV domain, the following ex- 
perimental conditions should be fulfilled: in case of the GRBs 
where keV-MeV photons are detected, the precision on At of 
the order of 10~ 4 to 10" 3 s is needed. For GRBs and AGN flare 
photons in the GeV-TeV range, the required precision should 
reach a level of a second for an individual photon. 

The procedures in use to determine time-lags rely on differ- 
ent statistical treatments of the data and most of them follow 
advanced procedures: 



• Cross Con-elation F unction (CCF) — BATS E (Band, 



1997b. H.E.S.S. dAharonian et all 120081) . Fermi 



(Abdo etal.[[2009albl) : 



Energ y Cost Function (ECF) — MAGIC (lAlbert et al. . 



2008), Fermi dAbdo et al.Lr2009allbh : 



Wavelet Transfo r ms — BATSE, HETE-2 , SWIFT 
fellis et all 120031 120061: iBolmont et all l2008h . H.E.S.S. 



(Aharonianet al. 



2008); 



Likelihood fit of QG model parame ters — INT EGRAL 



dLamon et all 120081). MAGIC (Martin ezand Errandol 



2009), H.E.S.S. IXbramows ki et alll201lh : 



Figure of merit formula (Eq. Q]) 
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Figure 2: Evolutio n of x 1 function of EL ,, for maxima and for minima found by 
CWT procedure in Bolmont et al. (2008). Each curve corresponds to a different 
choice of the energy bands, and then to a value of AE. All scenarios fulfill the 
condition ^(E^/ndf < 2. 



Figure 3: 95% CL conto urs for a and b f rom the 2-parameter fits for the fourteen 
scenarios considered by Bolmont et al. 12008). The box at the bottom left of the 
plots shows the position of contour centers. 



For robust results, the use of at least two methods which 
probe different aspects of the light curves is recommended. In 
addition, careful error calibration studies by Monte Carlo sim- 
ulations are mandatory for the extraction of the limits. 

3. Present data and results 

In this section the studies performed in view of search for 
Lorentz Symmetry breaking with Pulsars, GRBs and AGNs re- 
spectively are presented, following the best methods in use in 
this field. The analyses performed by the High Energy Tran- 
sient Explorer (HETE-2), by Fermi GBM/LAT mission and by 
the H.E.S.S. and MAGIC experiments are good examples of 
current studies in presently running experiments. The quality 
and the significance of the results will be related to the acquired 
statistics and the type of method in use. 

3.1. Pulsars 

The pulsed emission from the Crab pulsar and nebula has 
been first studied by the EGRET experiment JKaaren, 1 19991) 
with respect to Quantum Gravity scale limit. Following to 
the alignment of the Crab pulsations (with 33.18 ms period) 
from radio to X-ray and gamma energies, a 95% CL limit of 
1.8 x 10 15 GeV was set on Eqg- This interesting result was 
achieved due to an excellent sub-millisecond precision of the 
EGRET clocks. A more stringent limit from the Crab pulsar is 
expected in future from the low energy threshold data acquired 
by the MAGIC collaboration. 

3.2. Gamma-Ray Bursts 

The GRBs are the most violent phenomena observed in the 
Universe, detected as sudden and unpredictable bursts of opti- 
cal, hard X and y-rays, lasting tenths of seconds and presenting 



variable and unclassified light curves. The GRBs are of cos- 
mological origin and result from the death of a massive star 
or from the collapse of compact binary objects. Their very high 
variability of the order of milliseconds in the large energy range 
and distances going up to redshift values of 8, place them as ex- 
cellent candidates for searches of the non-standard effects in the 
photon propagation. 

An example of a typical analysis in the field in keV-MeV 
energy range has been performed by HETE-2 experiment. This 
mission was devoted to the study of GRBs using soft, medium 
X-rays and gamma-rays with instruments mounted on a com- 
pact spacecraft. The analysis of the 15 GRBs with measured 
redshifts collected by HETE-2 mission in ye ars 2001- 2006 



follow ed the procedure described in details by iBolmont et al 
(2008) and references therein. After the determination of the 



GRB interval time describing the start and the end of the burst, a 
cut above the background delimited the signal region to be stud- 
ied in further analysis. The originality of the performed analy- 
sis was the choice of various energy ranges where the time-lags 
were computed with tagged photon data provided by the FRE- 
GATE sub-detector of the HETE-2 for each GRB. The light 
curves of the 15 GRBs are shown in Fig.fl] It appears clearly 
that the signal-to-noise ratio decreases for large redshifts. 



Following Ellis et al. (2003), after a de-noising procedure 
of the light curves by a Disc rete Wavelet Transform (DWT) 
(Don oho and Johnstond.ll994t) . a search for the rapid variations 
(spikes) in the light curves for all energy bands using a Continu- 
ous Wavelet Transform (CWT) dMallafl, 1 19991) was performed. 
As a result, a list of minima and maxima candidates was fol- 
lowed by their association in pairs. 

Fig. [2] presents the evolution of;^ 2 (EQG) as computed with 
equation [10l around its minimum x 2 ■ (EoG)/ndf for maxima 
and mini ma together. The two -parameter linear fits, as dis- 



cussed by IBolmont et alj (|2008), show a somewhat different 
behavior in case of the maxima and the minima. However, no 
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Figure 4: Time-lags from high intensity GRBs as collected by BATSE, HETE-2 
and SWIFT with rescaled errors, which lead to robust results, EL~ > 2.1 X 10 16 
GeV at 95% CL fcllis et all[2006ll2008h . K is defined as in Eq.[9] but without 
the factor (1 + z) in the integral. 



significant preference of any value of Eqg is observed for most 
of the scenarios when considering both types of extrema. 

Concerning possible source effects, it has been known for a 
long time, that the peaks of the emission in G RB light curves 



are shorter and arriv e earlier at higher energies dFenimore et al 



1995; Norris, 2002). These intrinsic lags, which have a sign op- 
posite to the sign expected from LIV with subluminal photons, 
have a broad dispersion of durations, which complicates the de- 
tection of the Lorentz Symmetry violation effects. In the study 
with HETE-2 GRBs, a universality of the intrinsic source time- 
lags has been assumed. The 95% CL contours for a and b from 
2-dimensional fit are presented in Fig. [3] showing that both pa- 
rameters representing the time-lags expected from the propaga- 
tion and the source effects are strongly correlated. In this figure, 
the values of the offset parameter b are compatible with zero for 
all energy scenarios. Both fit results suggest no variation above 
2cr, so that the 95% CL lower limits of the order of 4 x 10 15 
GeV on the Quantum Gravity scale were derived. 

In conclusion of the "low energy" studies, the time lags 
from GRB emission detected by HETE-2, BATSE, RHESSI, 
SWIFT and INTEGRAL space-borne experiments, as shown 
in Fig. [4] for the measured time-lags, in the energy range of 
keV-MeV, provide robust lower limits on the Quantum Gravity 
scale of the order of 10 16 GeV as analysed in (lEllis et all 12006. 
2008h . 

More recently, the Fermi collaboration (lAbdo et al. 



2009a.b) has published two results on the Quantum Gravity 
scale lower limits reaching or exceeding the Planck scale from 
data taken with GBM and LAT sub-detectors, when detecting 
a powerful signal from GRB 0809 16C with a redshift value ~4 
and from a short burst (GRB 090510) at redshift of ~ 0.9. The 
light curves as measured in energy bins are shown in Figs. [5] 
and [6] In both cases, substantial time delays of several seconds 
were measured between photons detected at keV (GBM) and 



GeV (LAT). As there is no strong evidence for a separate high 
energy component, the entire energy range was used to derive 
the limits on the LIV linear term of E' > 1.5 x 10 18 GeV for 
the first GRB and E' > 1.5 x 10 19 GeV at 95% CL for the 
second one, attributing the measured lags to the source emis- 
sion effects. These results constitute a breakthrough point in 
constraining the linear term models predicting Quantum Grav- 
ity scale below Planck scale. However, it should be noticed that 
both results rely on individual photons with energy > 10 GeV 
which determines the energy ranges. The limits were computed 
either by simply dividing A(t) by the distance formula, or by 
CCF and Cost Function, with no statistical calibration proce- 
dure implied in the error estimation. To further progress in this 
domain of energy, future detection of several GRBs by Fermi 
mission with known redshift and measurable time-lags will be 
studied. 



3.3. Active Galaxies 

Blazars are variable AGNs, extragalactic sources producing 
y-rays via the gravitational potential energy release of matter 
from an accretion disk surrounding a Super Massive Black Hole 
(SMBH). Beamed emission, large inferred luminosities, rela- 
tivistic plasma jets pointing to the observer, and the flux varia- 
tions by large factors (flares) on hour scale in time, make them 
excellent objects for variability studies from radio to VHE pho- 
tons. In addition, the blazars are valuable transient candidates 
for searches of effects due to Lorentz Symmetry violation at 
Quantum Gravity scale. So far, more than 20 AGNs have been 
detected in VHE range with redshifts varying between 0.002 
and 0.4. In spite of the fact that several AGNs have been already 
detected, only few provide enough luminosity for the analyses 
of their variability and photon propagation effects. Here, the 
case of the observation of an exceptional flare of the PKS 2155- 
304 (z = 0.1 16) by H.E.S.S. is described in m ore detail. 

In 2006, the H.E.S.S. experiment detected (lAharonian et all 

20071) an exceptional flare of this source, with a high flux 
(10,000 photons recorded in 1.5 hours) and a high variability 
(rise and fall times of ~ 200 s during the night of the 28th of 
July). The over-sampled light curve of the flare is shown in 
Fig.|7]in two different energy bands. 

To measure the time lags between photons in two different 
light curves, two independent analyses were carried out, using 
two different methods: 



determining the position of the maximu m of the Mod i- 
fied Cross Correlation Function (MCCF, iLi et all (12004 ) 
which gives directly the value of the time lag. This method 
was applied to the over-sampled light curves of Fig. [7] in 
the energy bands 0.2-0.8 TeV and > 0.8 TeV, which cor- 
responds to a AE of 1 .02 TeV. 



using a Con tinuous Wavelet Transform (Mallafl Il999 : 



Bacrvl l2004h to locate with great precision the spike po 



sitions (extrema) in the light curves. An extremum of the 
low energy band was associated with an extremum in the 
high energy band to form a pair. 



o 
o 



o 
o 



£1 

"Si 



o 
o 



c 

1 



o 
o 



o 
o 




3000 



2000 ~£ 



1000 



1000 



o 
o 



500 m 



o 
o 



600 



400 ■» 



o 
o 



200 



o 
o 



o 
o 



40 60 

Time since trigger (s) 



Figure 5: Light curves in dif ferent energy bands for GRB 0809 16C as recorded by Fermi instruments GBM (the two panels at the top) and LAT (the three panels at 
the bottom). Plot taken from lAbdo et alj<20093) . 



10- : 



1Q 3 



10 z 



10 




' • » • • :. 

ftT* ."".I.M,'.. „?...: •.*..♦. 






0.5 1 1.5 

Time since GBM trigger (1 a May 2009, 00:22:59.97 ut) (s) 



Figure 6: Light curves in different energy bands for GRB 090510 dAbdo et aill2009bT) . The top panel (a) shows the energy of each photon as a function of time. 
The vertical lines correspond to different hypotheses about when the highest energy photon (3 1 GeV) was emitted. These different hypotheses lead to limits on Eq G 
varying from 1.19Ep to more than 100Ep. 
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Figure 8: Energy Cost Function (ECF) as a function of the time delay. The 
maximum value of the ECF provides the measured time-lag, here for the linear 
model. Figure taken from Albe rt et alj J2008I) . 
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Figure 9: The time-lag measured as a function of redshift for three active galax- 
ies (Mkn 501, MKn 421 and PKS 2155-304). The results of the Whipple and 
MAGIC experiments are compared with the present H.E.S.S. results. 



The fit of the MCCF curve with a Gaussian plus a first de- 
gree polynomial allowed for the determination of the time-lag 
between low and high energy photons of 20 s. In order to eval- 
uate the uncertainties on this result, a detailed simulation of 
10,000 light curves was performed in each energy band varying 
the flux within the error bars, leading to a calibrated error of 
28 s. The obtained time-lag being found to be compatible with 
zero, a 95% CL upper limit on the linear dispersion of 73 s/TeV 
was set. 

With the CWT method, two pairs of extrema were obtained 
giving a mean time delay of 27 seconds. A similar method as 
the one used for the CCF was used to determine the errors and 
they were found to be in a range between 30 and 36 seconds. 
A 95% confidence limit of 100 s/TeV was obtained for the lin- 
ear correction to the dispersion relations. The MCCF method 
leads to a limit of the E' > 7.2 x 10 17 GeV and the CWT to a 

confirmation of this value with a limit of the EL, > 5.2 x 10 17 

go 

GeV. The work in progress with Likelihood fit of the time-lag 
with individual phot ons provides a prel i minar y lower limit of 



E' QC > 2 x 10 18 GeV (lAbramowski et all 1201 ll) . thus approach- 
ing considerably the limits found with Fermi GRBs. 

The presented analysis with PKS 2155-304 flare can be 
compared with those performed by Whipple and MAGIC ex- 
periments with other AGN flares: Mkn 421 and Mkn 501. The 
Whipple collaboration has set a limit of 4 x 10 16 G eV using a 
flare of Mkn 421 (z = 0.031) feiller et all Il999h . More re- 



cently, the MAGIC collaboration obtained a lim it of 3.2 x 10 n 
GeV with a flare of Mkn 501 at z = 0.034 dAlbert et al.ll2008l) . 
Fig. Upresents the ECF plot as obtained bv l Albert et al. (2008) 
assuming a linear dependence on the Quantum Gravity scale. 
A similar value was derived from the same data with the likeli- 
hood fit method ([Martinez and Errandoil2009l) . Although a nice 
minimum has been obtained in the likelihood function, a deci- 



sive conclusion could have not been established because of a 
limited statistics collected with Mkn 501 flare. 

The results obtained with H.E.S.S. are more constraining 
due to the fact that PKS 2155-304 is almost four times more 
distant than Mkn 421 and Mkn 501, and the acquired statistics 
were higher by a factor of ten. Fig. [9] summarizes the results 
obtained so far with AGNs. 

4. Summary and discussion 

Quantum Gravity phenomenology applied to studies with 
astrophysical sources has known a growing interest in the past 
decade, especially since it was argued that Lorentz Symmetry 
could be violated or at least distorted by the quantum nature 
of space-time, with measurable effects on photon propagation 
over large distances leading to a modification of the light group 
velocity. 

The y-ray dispersion studies in vacuum with GRBs and 
AGNs provide the cleanest way to look for Planck-scale modi- 
fications of the dispersion relations at the first order level. The 
results obtained with these sources for linear and quadratic cor- 
rections to the dispersion relations (translated into limits on pa- 
rameters EL G and Eq G ) are summarized in Table 1. 

As discussed in this paper, the GRBs are good candidates 
for time of flight studies. The GRBs are easily detected by satel- 
lite experiments at very high redshifts up to z ~ 8 and up to a 
few hundreds of GeV in AE. Population studies have alr eady 
been carried out ( Bolmont et all l2008t Ellis et all [2006) and 
lead to limits of the order of 10 16 GeV. Th e best and most spec - 



2009b). 



tacular limit so far is E' G > 2 x 10 19 GeV (lAbdo et al 
The measurements of the AGN flares comes in complementar- 
ity to the GRB studies and yield at present lower limits on the 
linear term dispersion parameters only a factor of few below the 
Planck scale. 




Figure 7: Light curve of the PKS2155-304 flare during the night of 28th July 2006. The range of the energy is 200-800 GeV for the top panel and > 800 GeV 
for the bottom panel. The original data points (in black) are binned in 2-min time intervals. Th e zero time point is set to MJD 53944.02. Gray points show the 
over-sampled light curve, for which the 2-min bins are shifted in units of five seconds. Taken from Aharonian et al. 12008). 



The existing synergy between studies done in Fermi exper- 
iment and ground based y-ray telescopes opens new era in the 
presented studies. In principle, a consistency of the linear and 
quadratic term results derived with data from sources of differ- 
ent nature would provide convincing results on LIV phenom- 
ena. The coherence of the results in principle could be postu- 
lated. Alternatively, apersistence of different results with GRBs 
and AGNs may signify energy dependence of the effect. Com- 
bined analysis of GRB and AGN data may also lead to a better 
understanding of the source induced part of measured time-lags 
as it concerns analyses of sources at different redshifts. It has 
to be underlined here that a combined analysis needs adequate 
analysis methods and procedures for sources with different vari- 
ability and different experimental setups. 

Concerning the interpretation of the results, the most con- 
servative one is to attribute the observed delays from MAGIC, 
H.E.S.S. and Fermi to the source emission effects having for 
its origin low energy gamma accelerated by electrons and high 
energy gamma by protons. The source delays are different 
for GRBs and AGNs (also different for short and long GRBs) 
and intrinsic fluctuations in the relative emission times for high 
and low components for sources in a given class or individ- 
ual source have to be taken into account and have to be con- 
trolled as this aspect limits the sensitivity of these probes of 
the LIV effects. The only way at present to minimize these 
effects are population studies of GRBs and AGNs as a func- 
tion of redshift which still need modelling of the source ef- 
fects. The other important point concerns a potentiality of 
these studies to detect a propagation induced time-lag: this 
would be a much more difficult and challenging task due to 
the complexity of the studied emissions. In particular, there 
is no indication from the theory if the deviations are expected 
to be of the sub-luminal or super-luminal type, thus combin- 
ing in different way with the source emission lags. Various 
proposals of studies are currently worked out as presented by 



Amelino-Camelia and Smolinl (12009b . 

In conclusion, till now, no significant result on Lorentz 
Symmetry breaking has been obtained, analyzing large energy 
range for photons emitted by GRBs or Very High Energy data 
from flares of AGNs. The presented studies have shown a good 
sensitivity to the linear term in the photon dispersion relations, 
setting the lowest limit for the Quantum Gravity scale around 
10 19 GeV. To constrain at the same level the quadratic term 
would require the energies in the range of the Ultra High En- 
ergy Cosmic Rays (UHECRs) or extra-galactic neutrinos. In 
future, the redshift dependences will be explored to distinguish 
between intrinsic effects to the source or induced by Lorentz 
Symmetry Breaking. Further observations of both a high num- 
ber of GRBs and of AGN flares will be necessary to give robust 
conclusions on possible propagation effects and possible devi- 
ations from Lorentz symmetry at high energies. Present and 
future experiments such as Fermi, CTAlj or AGISQ will greatly 
improve our capabilities in this area. 
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Table 1: A selection of limits obtained with various instruments and methods for GRBs AGNs, and the Crab pulsar. Limits obtained for linear (Eq G ) and quadratic 



(E*) corrections are given. 



Source(s) 



Experiment 



Method 



Results" 



Reference 



GRB 021206 
GRB 0809 16C 



GRB 090510 



9 GRBs 



RHESSI 
Fermi GBM - 



Fermi GBM + LAT 



BATSE + OSSE 



Fit + mean arrival time in a spike 

associating a 13 GeV photon with the 

trigger time 

associating a 3 1 GeV photon with the 

stall of any observed emission 

wavelets 



15 GRBs 


HETE-2 


wavelets 


17 GRBs 


INTEGRAL 


likelihood 


35 GRBs 


BATSE + HETE-2 + Swift 


wavelets 


Mrk 421 


Whipple 


likelihood 


Mrk501 


MAGIC 


ECF 
likelihood 



PKS 2155-304 H.E.S.S. 



MCCF 



wavelets 
likelihood 



■-(.a , 



> 1.8xl0 17 GeV 

> 1.3xl0 18 GeV 
>0.8xl0 10 GeV 

> LSxlO^GeV 
>3.()xl0 10 GeV 
>0.7xl0 16 GeV 

> 2.9x10 s GeV 
>0.4xl0 16 GeV 

> 3.2x10" GeV 

> 1.4xl0 16 GeV 

> 0.4x10" GeV 
>0.2xl0 18 GeV 
>2.6xl0 10 GeV 

GeV 
>5.7xl0 10 GeV 

GeV 
>0.1xl0 10 GeV 

GeV 
>2.1xl0 18 GeV 



E' Qa > 0.3x10 



nl 



Crab pulsar 



EGRET 



shift of pulsation maxima in diiferent 
energy bands 



E^ r >6.4xl0 10 GeV 

QG 

E' >0.2xl0 16 GeV 



Boggs et al. (2004) 
Abdo et al. (2009a) 

Abdo et al. (2009b) 

Ellis etal. (2003) 

Bolmo nt et al. (2008) 
Lamon et al. (2008) 
Ellis et al. (2006. 2008) 
Biller etal. (1999) 
Albert et al. (2008) 

Martinez and Errando (2009) 

^haroniane^al. (2008) 

Abramowski et al. (201 1) 
Kaaret(1999) 



Note 
Tc 



/ 



a Results can be expressed as by Kostelecky and Mewes (2008) using the formulas F/ = 1/E QG and Y q = 3/2 ■ 1/(E QG ) 2 . 

'' Limit obtained not taking into account the factor (1 + z) in the integral of Eq.[3landl4l 

' The pseudo-redshift estimator i Pelangeon et al., 2006) was used. This estimator can be wrong by a factor of 2. 

d Only the most conservative limit is given here. 

e Photon tagged data was used. 

* The pseudo-redshift estimator i Pelangeon et al.. 2006) was used for 6 GRB out of 1 1 . 

g For HETE-2, fixed energy bands were used. 

; ' The limits of Ellis et al. 1 2006) were updated in Ellis et al. 1 2008) taking into account the factor (1 + z) in the integral of Eq. I3landl4l Only the limit obtained for a linear correction is 

given. 

' A likelihood procedure was used, but not on an event-by-event basis. 
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